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ABSTRACT 

Aims. We study the observational signature of flux emergence in the photosphere using synthetic 
data from a 3D MHD simulation of the emergence of a twisted flux tube. 

Methods. Several stages in the emergence process are considered. At every stage we compute 
synthetic Stokes spectra of the two iron lines Fe I 6301.5 A and Fe I 6302.5 A and degrade the 
data to the spatial and spectral resolution of Hinode 's SOT/SP. Then, following observational 
practice, we apply Milne-Eddington-type inversions to the synthetic spectra in order to retrieve 
various atmospheric parameters and compare the results with recent Hinode observations. 
Results. During the emergence sequence, the spectral lines sample different parts of the rising 
flux tube, revealing its twisted structure. The horizontal component of the magnetic field retrieved 
from the simulations is close to the observed values. The flattening of the flux tube in the photo- 
sphere is caused by radiative cooling, which slows down the ascent of the tube to the upper solar 
atmosphere. Consistent with the observations, the rising magnetized plasma produces a blue shift 
of the spectral lines during a large part of the emergence sequence. 
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1. Introduction 



Observational studies of magnetic flux emergence in the photosphere with the help of Stokes 
parameters started with the work of 



Brants 



and 



Zwaanetal 



( 19851) . Subsequent stud- 



ies showed that active regions grow by the successive emergence of magne tic flux fragments , 



1996; 



each containing only a small fraction of the total flux of the active region (IStrous et al 
Strous & Zwaarll999). Emerging fl ux regions (EFRs) of various sizes lead to magnetic co mplexes 



in a r ange of scales 



jHarveyl l 



19931) . from the largest active reg ions (up to 4 x 10 



1987b to small ephemeral regions (with fluxes down to 10 19 Mx; lHagenaar et all 199% . 



Mx; 



Zwaan 



While the initial studie s were restricte d to spectropola rimetric observations of Stokes-/, Stokes 



V and continuum images (IBrants 



1985a.b: 



a 



Strous 



19941) . the full Stokes vector has been used by 
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Lites et al 



Lagg et al 



(1998b and 
d2004 



Kubo et al 



(12003k for photospheric studies and by 



Solanki et al 



d2003l) and 



2007) for chromospheric investigations. Recently, observa tional studies of flux 



emergence with high spectral and spatial resolution have been carried out (iCenteno et al 



Cheung et al 



2008 



Ishikawa et al 



2008 



Martinez Gonzalez et al 



2007; 



Okamoto et al 



2007 



2008) 



with the spectropolarimeter (SP: lLites et al.l200lb whic h is part of the focal-plane package of the 50 
cm Solar Optical Telescope (SOTi rTsuneta et al]2008 ) onboard the Hinode spacecraft. The emerg- 
ing flux regions show horizontally oriented flux elements, which are generally regarded as the tops 
of emerging loops. The reported field stren gth of these ho rizontal flux elements varies somewhat 
among observa tional studies : 500 ± 300 G (Brantslll985alb ): 20 < \B\ < 600 G ( Lites etall ll998 



Sigwarth 



2000); 400 < \B\ < 700 G (Kubo et al. 2003); 650 G (Ok amoto etal 



(Brantsll985alb 


: Strous etal.il 996): 


km/s ( 


Okamoto et al. 


2008 


). 


Lites (2005) 


and 


Okamoto et al. 



20081) . The rising 
city are : 0.5 km/s 

al.l ll996l) : ~ 1 km/s dLites et alJU998h : < 1 km/s dKubo et alJboO^ : 0.3 



(2008) h ave shown t hat so me flux elements emerging in the 



Wiegelmann et al.l (12005k. who fou nd newly emerged loops 



Okamoto et al 



(2008 



2009) (hereafter OK2008 and 



photosphere exhibit a twisted field (cf . 
to harbour considerable electrical currents). 
OK2009) reported observations of the emergence of a coherent helical flux tube. 

A number of theoretical studies has addressed magnetic flux emergence from the convec- 
tion zone to the solar atmosphere. While some approaches treat the convection zone and over- 
lying atmosphere as plane-parallel hydro static layers (e.g. 



2003 



Archontis et al 



Cheung et alJl2007 



2004 



MagaraT2006l etc), others include the effects of convective flows (e.g. 



Shibata et al 



1989; 



Abbett & Fisher 



Martinez-Svkora et al. 2008 



Isobe et al.! 2008) 



In this paper we present for the first time detailed synthetic observational properties from 3D 
MHD simulations of an emerging, initially horizontal, twisted flux tube including its interaction 
with granular convective flows. We perform spectral line synthesis and inversion at spectral and 
spatial conditions similar to Hinode's SP/SOT. The results are discussed and some emergence prop- 
erties are compared with the observations of OK2008 and OK2009, even though it is clear that there 
are limits to such a comparison, arising from the fact that observations and simulations do not occur 
in the same solar environment. Nonetheless, we believe that the similarities or differences revealed 
by a comparison will provide fresh insight into both simulations and observations. 



2. Flux emergence simulations 



2005) 



The simulation discussed here has been carried out with the MURaM code (IVogler et al 
which solves the compressible MHD equations taking into acount non-local energy exchange via 
non-grey radiative transfer as well as the effect o f partial ionization on the equation of state. Details 
of the simulation setup have been described by lCheung et al.l (12007b . The top and bottom bound- 
aries are open and allow for mass transfer. The magnetic field is matched to a potential field at the 
top boundary. The initial condition consists of a horizontal buoyant magnetic flux rope (i.e., a flux 
tube with a twisted magnetic field) embedded in the convection zone at a depth of 1.35 Mm below 
the solar surface. The initial axial field strength and characteristic radius of the tube are #0 = 8500 
G and Ro = 200 km, respectively, which gives a longitudinal flux of O = 10 19 Mx. The dimen- 
sionless twist parameter of the tube (corresponding to the ratio of the transverse and longitudinal 
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Fig. 1. Maps of the magnetic field vector retrieved from Stokes inversions. The colour code indi- 
cates the amplitude of the vertical (line-of-sight) component of the magnetic field. The orientation 
of the horizontal field components is shown by arrows. The panels a to d correspond respectively 
to: t = 11.3, 13.0, 16.1 and 23.4 min. The blue contour in the upper-left panel outlines regions 
where the horizontal component of the magnetic field is larger than 300 G. The vertical black line 
in the upper-left panel indicates the location where vertical cuts in the MHD cubes are plotted in 
Figs.[2]and[3] 



fie ld strengths at a dista nce r = Rq from the tube axis) is A = 0.5. This setup corresponds to run Ul 



of 



Cheung et al.l (12007k . To ensure that the temperature structure in the photosphere is sufficiently 



realistic for the synthesis of Sto kes parameter s, we have repeated this run with non-grey radiative 
transfer using fou r opacity bins ( Vo^leJl2004 ). 

As shown by ICheung et al.l (12007b . the setup described above allows a magnetic tube with a 
flux of 10 19 Mx to overcome the downward directed drag of convective downflows and emerge 
as a coherent structure, although the effects of the convective pattern are still impressed upon the 
emerging flux pattern at the solar surface. 



3. Spectral line synthesis and Milne-Eddington inversion 

We have computed the Stokes parameters for the iron lines Fe I 6301.5 A and Fe I 6302.5 A 



This has been done in LTE using the STOPRO code in the SPINOR package (iSolanki 



1987 



Frutiger et al.ll200Qh . For the atomic para meters characterising these two lines a nd the iron abun 



dance we have used the same values as 



d2005aL : 



Borrero et al 



(20061) and 



Shchuki na & Truiillo Buenol (1200 II) : 



Shelvag et al] d2007l) . 



Khomenko et al 



The resulting Stokes data are degraded to Hinode's SP/SOT resolution employing an optical 
model of the tel e scope and taking into account a defocus value typical of the Hinode's SP following 



Danilovic et al 



(2008). We use a pixel size of 0.32 arcsec. In line with typical characteristics of 
Hinode SP/SOT, we degrade the Stokes spectra to a spectral resolution of 30 mA and a spectral 
sampling of 21.6 mA. A photon noise of 10" 3 I c is also added. 
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The degraded synthetic Stokes data are inverted using a Milne-Eddington inversion procedure. 
This is done by using the Helix code (ILagg et al.ll2004l) . Both spectral lines are inverted simul- 
taneously. The existence of two atmospheric components within each pixel was allowed for, one 
harbouring a magnetic field, the other being field-free. The fraction of the surface area occupied 
by a magnetic field within each pixel is the "filling factor". The inversion procedure searches for 
the optimal combination of parameters corresponding to the minimum of the merit function. The 
following free parameters are considered: the magnetic field strength, its inclination and azimuth, 
the ratio of line center to continuum opacity, the Doppler width, the damping constant, the slope 
and surface value of the source function, the Doppler shift (giving the line-of-sight velocity) and 
the filling factor. 



4. Results and discussion 

Figure [T] shows a time sequence of the magnetic field maps retrieved by the inversion. Clearly, the 
emerging field expands laterally with time. Similar expansion has been reported in observations of 
the emergence of a flux rope in an active region by OK2008 and OK2009 (see Fig. 2 in OK2008 
during the first half of the rise of the observed flux rope). 

The orientation of the horizontal components of the field is indicated by arrows. The twist of the 
emerging flux tube is revealed by the orientations of the field lines, although the apparent structure 
is deformed by the influence of the granulation. This orientation changes along the emergence 
sequence (see the difference between the two snapshots at t - 1 1.3 min and t = 16.1 min). 

At t = 11.3 min the spectral lines sample the upper part of the emerging flux tube, where the 
field lines are predominantly oriented from upper left to lower right in the plotted frame. Later, at 
t - 16.1 min, the flux tube has emerged further and the spectral lines sample its lower part, where 
the field lines in this twisted structure are oriented from the lower left to the upper right. A change 
in the field orientation has been reported in the observed event by OK2008 and OK2009 although 
the observed event takes place in an active region. 

The mean value of the horizontal component of the magnetic field in the region outlined with 
the blue contour (near the polarity inversion line) in the upper left panel of Fig.[T]is 550 G at t = 11.3 
min. This value is close to the one (650 G) reported by OK2008 and OK2009. The horizontal field 
strength exhibits higher values in the close vicinity of the polarity inversion line and decreases with 
time along the emergence sequence shown in Fig.Q] The upward velocity of the emerging flux rope 
in the same region outlined in blue is 670 m/s. This value is larger than the 300 m/s reported in 
OK2008 and OK2009 and indicates a faster emergence rate of the simulated flux rope. The mean 
value of the magnetic field inclination with respect to the vertical in the region outlined in blue is 
nearly zero degrees (zero corresponds to a horizontal field) with a standard deviation of about 30 
degrees. In OK2009, it is reported that the mean inclination near the polarity inversion region is 
nearly zero, with a scattering of the order of ±10 degrees. The mean vertical magnetic field in the 
same outlined region is nearly zero. The one reported in OK2009 is ±200 G. 

The total unsigned magnetic flux through the maps in Fig. [T] is calculated from the inversion 
results as: Flux = A ^ \B z j\fi, where A = pixel area, B z i = B z in pixel i and = filling factor 

i 

in pixel i. For the four snapshots a-d in Fig. [T] we obtain the following values: 1.15 x 10 20 Mx, 
1.42 x 10 20 Mx, 1.45 x 10 20 Mx and 1.05 x 10 20 Mx, respectively. The increase of the unsigned 



Yelles Chaouche et al.: Flux emergence simulation and comparison with observations 5 

flux in the first three snapshots is due to flux emergence. Subsequently, the decrease is due to flux 
cancellation of opposite-polarity flux elem ents. These element s are formed by the fragmentation of 
the tube under the action of flux expulsion ( Vogler et al.ll2005 ). The encounter of opposite polarity 
flux elements in the intergranular lanes leads to the decrease of the unsigned flux. 

2D vertical cuts through the MHD cubes perpendicular to the axis of the original flux rope are 
shown in Figs. [2] and [3] at t= 11.3 min and 16.1 min, respectively. They lie along the horizontal 
coordinates represented by the black line in the upper left panel of Fig.Q] 

The rising flux tube flattens considerably after reaching the solar surface. The upper parts of the 
flux tube (above the solar surface) rapidly cool down radiatively. Therefore, the plasma becomes 
denser and loses its buoyancy. On the other hand the bulk of the flux tube keeps rising. This is 
due to the fact tha t the rising sub- surf ace parts of the flux tube keep pushing the plasma upward 
above the surface (iCheung et al.ll2007l) . The material above the surface also expands laterally due 
to the strong stratification of the photospheric layers and partly descends back into the convection 
zone at the sides of the flux tube. We see a blue shift during most of the emergence sequence at the 
height where Fe I 6301.5 A and Fe I 6302.5 A are formed, except for downflo w lanes along the 



flux rope where we find a red shift (ICheung et al 



2007 



Yelles Chaouche 



2008). A blue shift has 



been observed as well by OK2008 during a large part of the emergence sequence. 

Figure [2] shows that the flux rope has a relatively regularly twisted core with a field strength 
approaching 1 kG near and below the solar surface. Note that the Stokes- V signal is formed over a 
height-range sampling mainly sub-kG field. Convective flows affect the outer parts of the emerging 
flux tube more strongly, so that the field configuration becomes more complex than in the central 
part. This is best seen in the field inclination (lowest panel of Fig. [2]) and is mainly restricted to 
the sub- surface layers. At the height of the line formation the field is rather homogeneous and the 
opposite polarities are separated (see also the top left panel in Fig.[T]). 

At t = 16.1 min (Fig. [3}, one can still see the downflows beside the flux tube, where the field 
is predominantly vertical. The flux tube has become wider, and there are patches of downflowing 
plasma above the solar surface in the body of the flux tube (e.g. near the abscissa x = 5500 km). 
The horizontally oriented core of the flux tube has greatly shrunk. Part of the flux has been advected 
to downflow lanes and concentrated to kilo-Gauss field strength. This produces a strong circularly 
polarized signal correlated with downflow lanes in the emergence area. At this point in time the 
flux rope is starting to loose its coherence even at the height of the line formation (Fig. [3}, with 
multiple islands of both polarities being present. 

By the time of the last plotted snapshot (Fig.QJl) many islands of opposite polarities are visible. 
Nonetheless, a predominance of positive polarities in the upper and negative polarities in the lower 
part of the frame can be discerned. In vertical cuts similar to Figs. [2] and [3] we show how the 
inhomogeneity of the magnetic structure slowly propagates upwards. 



5. Discussion 

We have studied the observational signature of a twisted flux tube emerging through the solar 
photosphere with the help of 3D MHD simulations. This emergence event reproduces some of the 
characteristics of the event studied by OK2008 and OK2009. 



6 
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Fig. 2. Vertical cut through the simulation box at the location shown by the black line in the upper 
left panel in Figure Q] at t = 11.3 min. The panels represent, from top to bottom, the temperature 
fluctuations with respect to the horizontal average, vertical velocity (blue represents upflow), mag- 
netic field strength and magnetic field inclination. The solid black lines indicate the location of 
T5000 = I- The dash-dotted lines show the height at which the contribution function of Stokes-/ 
reaches its maximum, and the dashed lines indicate the same for Stokes- V. 



There are also major differences, however, which we now briefly discuss: Whereas the channel 
of emerging flux closes on itself according to OK2008 we find no sign of this in our simulations. 
Instead, the emerged field at the surface continues to expand. At the same time the structure of 
the vertical component of the field becomes i ncreasingly inhomog eneous. This is similar to the 
evolution of an emerging loop system studied by lCheung et al.l (12008b . Also, the simulated flux rope 
emerges within a few tens of minutes whereas the observed emergence takes place over roughly a 
day. 



Cheung et al.1 (120071 120081) have made a parametric study where they discuss the dependance of 



the emergence properties on the initial magnetic field strength and twist. The speed of emergence 
depends on the total buoyancy of the tube (related to total flux and total density deficit). In that 
sense, more buoyant tubes rise more quickly and can modify the downflow pattern. Very buoyant 
and twisted tubes tend to modify the granulation pattern in an obvious way that is lacking in the 
OK2008 event. We can then assume that the buoyancy of the observed emerging flux is not very 
large, causing no anomalous granulation and has gentle flows. This is consistent with the fact that 
the observed emergence has an upflow speed lower than the simulated case and takes place over a 
longer period. We know the behaviour of the simulated emergence onl y over a rather limited rang e 



of parameter space (initial magnetic field strength and twist) studied by lCheung et al 



(12007 



2008), 



which is still different from the regime of the OK2008 case. In future work, we plan to investigate 
flux emergence leading to and within active regions. The emergence process inside a pre-existing 
field will be of particular interest. 
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Fig. 3. Similar to Figure[2] at t = 16.1 min. 



